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ABSTHACT 


In tvo preceding papers (Wu and Lee, 1979; Lee et al . , 

1980) an effective loss-cone distribution function is used in the 
stability analysis relevant to the theory of auroral kiloaetric 
radiation. Furthermore, in these discussions it is assumed that 
the population of low-energy background electrons prevails over 
that of the energetic electrons. Obviously, this assumption may 
not be justifiable at high altitude. The present work con- 
siders a loss-cone distribution close to that observed and, more- 
over, inclu-les the effect of a parallel electric field along the 
auroral magnetic field lines. It is found that this effect tends 
to enhance the growth rate. We have also removed the assumption 
concerning the high population of background electrons. We find 
that at high altitudes, where the density of energetic electrons 
may dominate and the usual cold-electron dispersion relation is 
invalid, the cutoff frequency of the I mode can be significantly 
affected by the energetic electrons and thereby the growth rate 
is indirectly altered. In the present discussion, the emission 
of the auroral kilometrlc radiation is studied over a range of alti- 
tudes. have considered three density models and two electric 
field models. Temporal growth rates and e-folding distances of 
spatial amplification are computed and discussed. 
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I. rmoDUCTioN 


la two preceding artlclM (Wu acd L««, 1979; Lm «c al . . 1980) a 
th«ory of th« aurorsl Icilotuetrlc radiation (AK20 Is proposod. Th« basic 
physical notion is that only a fraction of the energet.Vc (VW) electrons 
originated in the plasma sheet can precipitate in the auroral region be- 
cause those with pitch angles greater than the atmospheric loss-cone angle 
are expected to be reflected by the convergent magnetic field before they 
con reach the upper atriosphere. Ihesc reflected electrons are expected 
to possess a loss-cone distribution function which can result in the 
amplification of electromagnetic waves with frequency very near the elec- 
tron gyro-frequency when the electron plasma frequency is sufficiently 
lower than the gyro- frequency (Wu end Lee» 1979). The emission mechanism 
is basically a cycle tron maser process. 

In addition to the hypothesis fhat the upgoing keV electrons can 
possess a loss-cone dlstribution.Wu end Lee (1979) also suggested that 
the parallel electric field of inverteo-V events (Gumetc, 1972; Mizera 
and Fennell, 1977; Croley et al . . 1978; Mozer et al . . 1977; Wescott ec al . . 
1976) can deplete the low-energy electrons of ionospheric origin in the 
auroral flux tube, thereby reducing the local value of the electron plasma 
frequency. This process can increase the gro^rth rate of the instability 
and, moreover, resu'.t in temporary wave trapping inside the source region. 
Evidently both consequences are Important to the emission process of AKR. 

Since the publication of the first article by Wu amd Lee (1979), we 
have been «.ncouraged by the facts chat Che S3-3 observations confirm our 
loss-cone hypothesis (private communications with Mlzera, Fennell and 
Croley) and the ISIS-1 results clearly show that the predicted density 
depletion in the AKR source region is real (Benson and Calvert, 1979; 
Benson, Calvert and Klumpar, 1980; Calvert. 1981). Furthermore, 
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oba«rv«elona by Bavkaya (Curnact and Craaa; 1978), by 7oyagar 1 
(Kaiaar at al . . 1978) and by ISIS>1 (Banaon and Calvaxt, 1979) all aaaa 
CO indlcaca that cha polarixatlon of AXR la In cha faac axtraordlnary noda 
(or X soda) .which la conalacant with our thaory. Our analyala ahowa chat 
although tha aapllflcaclon of cha ordinary iwda (or 0 soda) la poaalbla, 
cha procaaa la much laaa affactlva chan chat for cha X noda. Of couraa, 

It ahould ba nantlonad chat cha X ooda polarlxacion haa alao baan pra- 
dlecad by ochar chaorlaa (Malroaa, 1976; Grabba, 1980). Phyalcally It 
can ba aaally undaratood frcm llnaarlzad 71at.)w thaory why Cha growth rata 
aaaoclacad vlch cha X noda la largar chan chat of cha 0 noda. Furchamora, 
cha X noda vavaa with fraquanclaa cloaa co cha cyclotron fraquancy nay 
hava group valoclclan much snail ar chan cha apaad of light, a faatura 
dlffaranc fron chat of 0 noda. Consaquantly, cha spatial aapllflcaclon 
of cha X noda la aora affactlva chan cha 0 noda, avan If both hava tha aame 
caaporal grovch raCas. 

A shorCconlng of cha pracadlng works (Uu and Laa, 1979; Laa at al . . 

1980) was polncad out by Vu (1980) . In ordar Co danonacrata chat 
an alactromagnaclc cyclotron Inawablllcy can axlsc va htva consldarad an 
affactlva losa-cona distribution which varlas with according co cha 

following functional form 

y *\i 

Expraaalon (1) Is fraquancly consldarad In plaana physics co Invaatlgata 
scablllty problama concamlng chargad parclclas conflnad In a nagnatlc 
nlrror. Tha Indax 2.(2. ) 1) uaasuraa cha staapnass of cha gradlanc 
of an affactlva losa-cona dlscrlbutlon. Such a dlacrlbution 
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li, hov«v«r» fllgnlflcancly diff«r«ac from obo«rv«d wtatcosphorlc loss- 
coas dlstrlbuclons vhlch srs dlscusssd In Ssctloa II. Ths qussclon which 
srlsss asCursJLIy is whsehsr chs two typss of dlstrlbuclon fvmctlons can 
load to qualieatlvsly dlffsrsnt conclualoas (within ths contsxt of llnair 
atabllity anslysaa) . Furtharmors, ths rsallstic loss-cons distribution 
can bs significantly aodifisd whsn thsrs is a parallal slsctrlc flald 
bslow or at chs sics of obssrvation. Ths slsctrlc flsld tsnds Co broadsn 
ths loss-cons nsar v„ ■ 0 and chsrsfors can affscc chs calculation of 
Chs tsmporal growth rats. 

Anoclsr gsnsrallzacion of chS prscsdlng works is dssirabls. That is, 
prsvlously vs havs assumsd chat chs dsnsicy of ths backscactsrsd slsccrons 
is much highsr chan chat of chs ksV slsctrons. This assumption snablsa 
us to caleulats chs dispsrslon rslatlon by ignoring chs snsrgsclc slsctrons, 
which only contrlbuts to ths amplification of Chs radiation. Howsvsr, at 
high altitudss ( ) 1 fsv backscactsrsd slsctrons (sxcspt choss Chat can 
attain snsrgiss in chs ksV rangs) can rsach that hsighc without bsing 
rsflsctsd by ths dscrsaslng slsctrostatic potential assoclatsd with chs 
parallel slsctrlc flsld. (Hers vs havs implicitly postulated chat ths 
potential drop below the alcicuds of one is of the order of one kilo- 
volt or more.) If this conjecture is true, the assumption that chs density 
of chs backscattersd electrons is high must be relaxed. In this case both 
the dispersion relation and ths stability of the radiation art determined 
by the energetic electrons. Hence, the stability analysis becomes more 
complicated chan that presented by Lee et al . (1980) . 

In view of the aforamsneioned considerations, a series of investigations 
have been carried out. The purpose of this paper is to discuss these studies 
and to present the relevant results. Tae organization of the paper is 
described as follows. In Vec. II we discuss the physical model on which 
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ch« •cabilicy snalyslt it ca?ri«d ouc. Th« foraulAClon of ch« aCabillty 
chtor/ Is r sported in St^c. TII. Finally, Sec. 17 prsasncs aunsrical results 
end discussion. 


A 



II. PHYSICAL MODEL 


It ia 5.pproprlACa to flrat diacuaa tha phyilcal modal to ba coo- 
fldarad bafora analyzing tha atabllity problam. Thara ara aavara.’ baaic 
conaidaratr.ona which ara diacuaaad aapar^italy aa followa. 

A. Modal Dlatributlon F »-i nctlon Baaad on Spacecraft Obaarvatlona 

Tha moat applicabla obaarvatlonal reaulta available now ara thoaa 
acquired with Che S3-3 satallita, which haa bean operative in tha apparent 
source region of AIQ. Particle dlatributlon functions ara measured by the 
S3-3 spacecraft in the auroral acceleration region at various altitudes 
below about 8,000 kilometers (Mlzera and Fennell, 1977; Croley et al . . 

1978; Fennel and Comey, 1980). Because the in situ measurements of the 
auroral electron dlatrlbution functions provide valuable information which 
can be the basis of our theoretical model, a brief discussion of the ob- 
setvationi'l results is not only relevant but also pertinent. 

Figure 1 shows typical invettad-V signatures as observed by S3-3. 

This particle-wave spectrogram from Day 237, 1976 shows electron and ion 
differencial energy flux from 0.2 - 30 keV and .09 - 3.9 keV, respectively . 
Also shown In Che upper panel is the AC electric field amplitude (Mozer 
•C al . . 1977) with frequenci-s of .CJ - 100 KHz. (Note chat S3-3 is not 
capable of making in situ obaervaclons of AKR because the high-frequency 
cutoff of the AC electric field detector lies below the local electron 
gyro frequency everywhere along the spacecraft trajectory.) In this example 
Che spacecraft encounters several invert ed-V structures between 69.8* and 
74.5* invariant latitude. Electron acceleration to several kilovolts and 
ion acceleration to 1 - 3 keV are observed throughout this region. 
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For the purpose of illustration we present a distribution function 
which was leasured at 17473 seconds UT(Fig. 21.) The + V| axis corre- 
sponds to downgoing electrons and the - Vj| axis to upcoming electrons 
observed at the spacecraft. There are several interesting features. First, 
overall the distribution function is fairly isotropic outside the loss-cone 
region. Second, the loss-cone region for the upcoming electrorj is parti- 
ally filled. Third, the loss-cone distribution appears to be modified by a 
parallel electric field. Fourth, a population of trapped electrons appears 
to be in exis ence. [The definition of trapped electrons may be found in the 
article by Chiu and Schulz (1978)]. 

ii- sed on the observations, v/e consider in the subsequent discussion 
a model distribution function of the following form 


F - A 
e 



- 0 


for 


2 

V, 


2 ^ 
Vh + * 



for 


2 

V, 


Vn + 


( 2 ) 


where A is a normalization constant, ^ and b are two positive numbers 
which are determined by the potential difference and the mirror effect. 


4 = 


2lel(4j^ - 4) 


m 


B 

and b = — ^ 1 


(3) 


e 

( and the electrostatic potential and magnetic field at the 

mirror point, respectively), in the context of the inverted- V structure 
itself, the quantity ^ represents the potential drop between the atmo- 
sphere and the point of observation. This can be estimai;td by measuring 
the ion beam energy on the enhancement in the upcoming loss-cone angle 
(see Mizera et al. , 1981). The step-function approximation is justifiable, 
if the velocity spread at the loss-cone boundary is small in coupartson 
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Fig. 2 . An Electron Veloclty^Spece Dletrlbutton Taken at 17473 Seconda 
UT on Day 237, 1976. The +v,, axle correaponda to dovr^iolng 
part Idea and the axla correaponda to upcoming partlclea. 
Th» 100-km loaa cone, corrected for parallel potential drop 
beneath the point of ohaervatlon, la ahaded. 
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with th« random valocity. If wa raquira 

d^v F - I 

J • 

ic can ba shown chac 



At this point a coupla of con&ants ara approprlata. Obviously, 
wa hava not Ircludad tha trappad alactrons. Of coursa, if tha parallal 
alactrlc flald only occurs balow tha slta of obsarvatlon, tha crapping 
procass naad not ba consldac^ anyway. Tha main raason is for simplicity, 
Anothar faaCura which wa hava not Cakan Into account In Eq. (2) Is that 
tha downgoing alactrons soma time can form a baam. This pciat shoaVd 
not posa any problam, bacausa In our analysis wavas ara considarad to 
hava a parallal componanc of tha wava vactor k pointing upward and 
Chasa wavas can only rasonata with tha upgolng alactrons. In short, tha 
functional form of tha distribution with + v„ la unimportant for tha 
itablllty problam of Intarasc to us. For tha sama raason, tha fcatura 
of a doubla loss-cona dlsplayad by Eq. (2) should not concarn us avan 
though wa ara incaraacad In tha raal situation In which only tha upgolng 
alactrons poasasa a loss-cona distribution. 

D. Magnatlc Flald and Elactron Dansity 

Tha medal distribution function dascribad by Eq. (2) dots dapand 
upon tha local valua of b (b = ** * “ntus a magnatlc~f laid 

modal la naadad, bafora wa can procaad furchar. Flgura 3 dascrlbas a 
modal considarad In tha prasant study. Tha magnatlc flald otrangth is 
axprassad In tanni, of tha alactron cyclotron fraquancy, (-1 aj B/m c) . 
which is plottad varsus altituda. In tha sama flgura savaral local loss- 
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R (Earth Radii) 


Plot Showing the Dependence of Electron Gyrofrequency and Loaa 
Cone Angle on Altitude Along Auroral Field Lines. 






we have considered B ^ O.U Gauss. 

tnax 

Gonctmins t‘ie eleccron densities we consider three models. In the fol- 
lowing we restrict our discussion to the altitudes ranging from 4000 km to 
8000 km. Of course, doing this does not necessarily mean that we are assuming 
that AKR cannot occur outside this region. In fact ,ISIS-1 observations in- 
dicate clearly that generation of AKR can take place below 4000 km (Benson 
and Calvert, 1979. Censon et al . , 1980; Calvert, 1981). rever, the general 
opinion is that the 4 x 10 - 8 x 10 km altitude range along the auroral field 

lines probably marks the most active source region of AKR. In Table 1, Che 

three density models are described. In Che same table we also list the values 

2 2 

of the local electron gyro- frequency and the corresponding ratios 

where u denotes the local eleccron plasma frequency. In Fig. 4 we plot 
P« 

Che local values of and u as functions of Che altitude. At a given 

e pe 

altitude, there are three values of u for different density models. We 

note here that u) is determined by the total electron density (i.e., 
pe 

n^ + n^ , n^ and n^ are Che density of the backscattered electrons and 
that of the energetic electrons, respectively). 

C. Parallel Electric Field 

The parallel electric ff ild which exists in the auroral acceleration 
region has been suggested and observed. However, the typical structure and 
distribution of such a field or fields still represent a controversial subject. 
A variety of possibilities have been suggested, as reviewed by Shawhan et al . 
(1978). Later the subject was further discussed by Chiu and Schulz (1978), 

Kan, et al . (1979), and many others from a theoretical point of view. 

Recent observations (c.f. Fennell and Gomey, 1980) suggest that the 
auroral potential drop can assume a variety of distributions in altitude. 

We therefore have considered some simple models. Since our lain interest 
is to investigate how the presence of a parallel electric field affects the 
generation of AKR, we have chosen to limit our study, for the purpose of this 
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Table 1. Density and Magnetic Field Models 


SfODEL 

ALTITUDE 

(km) 

“b 

(ca 

( CJB ) 

(kH*) 


b 5 J|i. 1 

a 8 


4,000 

40 

5 

373 

0.026 

2.5 


5.000 

25 

5 

281 

0.031 

3.7 

01 

6,000 

20 

5 

216 

0.033 

5.1 


7,000 

15 

5 

170 

0.056 

6.8 


8,000 

10 

5 

135 

0.066 

8.7 


4,000 

30 

5 

373 

0.020 

2.5 


5,000 

15 

5 

281 

0,020 

3.7 

02 

6,000 

10 

5 

216 

0.026 

5.1 


7.000 

5 

5 

170 

0.028 

6.8 


8,000 

3 

5 

135 

0.035 

8.7 


4,000 

4C 

5 

373 

0.003 

2.5 


5,000 

20 

5 

281 

0.005 

3.7 

03 

6,000 

0 

5 

216 

0.009 

5.1 


7,000 

0 

5 

170 

0.014 

6.8 


8,000 

0 

5 

135 

0.022 

3.7 ! 
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Fig. 4a. Plot Showing the Modeled Dependence of the Electron Plasma 
Frequency for the Altitude Range of Interest. 
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Altitude (km) 


Fig. 4b. SaiM as Fig. 4a, but Shoving the Dependence of the Ratio of 
the Squares of the Plasma Frequency to the Electron 
Gyrofrequency. 




paper, to two fairly extreme models. These two models are displayed in 
Fig. 5 and numerical values are given in Table 2. The first model assumes 
that the parallel electric field is distributed over a broad altitude range, 
while the second model assumes a narrow potential structure with potential 
drop only above 4000 km. 
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R (Earth Radii) 


Fig. 5. Plot of the TWO Potential Distribution Models used In the Wave 
Growth Calculations. Model 1 deplete a fairly broad potential 
distribution In altitude, while Model 2 contains a potential 
drop confined to ~3000 ka. 
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HI . D CSPERSION EQUATION 


In ord«T CO study ch« aaplificacion of chs rsdisclon w« wane co 
Invaatlgaca cha disparsion tqiiaclon which can ba darivad from llnaarlxad 
Vlasov aquaclona and la wall known in cha licaracura of plasoia physics. 

Tha ganaral form of cha disparsion aquaclon Is complicated if cha plasma 
is magnaclzad and Cha parcurbatlon has mixad alacCrcscacic and alecCro- 
magnaclc polarization. Furthermora, If tha wave fraquancy of incarast 
is closa co cha alaccron cyclotron fraquancy or Ics harmonics, tha rala~ 
tiviscic efface on cha mass of cha resonant alaccrons may be important, 
even for VeaV alaccrons. In this case, the solution of cha dispersion 
aquation is noc easily obtainable. The problem is sotLawhaC simplified 
whan Cha low-energy background electrons have a density do min a n t over 
that of the energetic electron#;. In this case tha real part of the wave 
frequency is pradetarminad by cha background electrons. Such an approxi- 
macion has bean considered by Lea at al. (1980), but it is not generally 
valid, as discussed earlier. Nevertheless, the approximate expressions 
for growth rates given either by Wu and Lea (1979) or by Lee et al . 

(1980), contain a piece of vary essential Information. That i.. , the 
resonant condition between the wave and electrons becomes qua'.itativ.4ly 
different when tha relativistic affect is included. A discussion of 
this point is not only relevant, but highly important. This is the 
purpose of Subsection III. A. 

Mathematically the relativistic effect also complicates the cal- 
culation of the dielectric censor. Usual.ly, when the relativistic effect 
is neglected, the integration with respect to v^ can be carried out 
for a certain class of non-equilibrium (unperturbed) distribution 
functions. The inclusion of the relativistic effect makes the 
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InccgraCion vary difficult. If aot lopoaslbl*. ^owavar, on* iiJiy adopt 
an approach dlffaranc from tha convantlonal nathod of calculating tha 
dlalactric tansor. Thla pola' la praaencad In Subnactlon 

Finally, tha dlaparalon aquation of intaraat to ua la praaancad 
In aubaactlon III. C. Tha aquation la axpraaaad in taraa of a nu^ar of 
intai^rala In T (tlaa) apaca. Evantually va hava to uaa ntimarlcal laathod 
to obtain tha aoluclon of thla dlaparalon aquation. 

A, Raaonanca Condition 

From Wu and Laa (1979) and Laa at al . (198C) , va taa that a raaonanca 
condition 


U) 

r 




N’',! 


- 0 


(5) 


playa an aaaantlal rola In tha datartJnatlon of tha tanporal growth raca. 
(Hara, tha approzlauttlon of vaakly ralatlvlatlc alactrona haa baan applied.) 
Thla condition can ba rawrlttan aa 


c 




N coa 


\ ^ w A A 

0 ) « -§• N coa 9 - 

/ 

a 



( 6 ) 


whara coa 9 ■ ki,/k ; N - kc/uj^ ti tha Indax rafracvlon; danotas 

tha wava fraquancy and • lajs^/m^c la tha alactroa cyclotron fra- 
quancy. Tha raaonant valocltlaa can ba raarllly plottad, aa ahown in Fig. 

6. Claarly, tha gaomatrical configuration la a clrcla which la cant ar ad 
at ■ c(u) N cot 9/n ) and haa a radlua R ■ f (tu^N^coa 9^/S^^) - 2(w /JJ -1)]^. 
Evldantly, tha dlacuaalon la maanlngful only whan 

2 

N'^coa 8 > 2 - ij- 

• ORIGJNAL . 

OF POOR QUALITY 
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-^N^cos‘e-2^ 


J!^ Cl) “ Cl) |t *" <^0 


Scheutttlc Plot Showing th* Shape of the Relatlvlatlc Resonant 
Surface for the Extraordinary Mode. The radius of the 
reaonitnce circle la ahown in the inset. 





Clearly, for u Q , the circle does not enclose the origin. One can 
’ re 

2 2 

easily envision that when N -ooand v « c , ttie resonant circle degene- 
rates into a straight line 

k cos e = - l) (7) 

which describes the usual Doppler effect and normal cyclotron resonance. 

The pertinent point here is that when u)^ , N cos6<l and v/c <1, 

all terms in Eq. (6) may be of the same order of magnitude. In that case 
the neglect of the relativistic effect in Eq. (5) is not justified. 

Moreover, it is seen that the inclusion of the relativistic term 
in the resonance condition tends to limit the resonant velocities to small or 


whereas the usual resonance condition described by Eq. (7) is independent 
of v^ , meaning electrons with arbitrary can interact with the wave 
as long as satisfies Eq. (7). Consequsntly, for a given unperturbed 

distribution the two resonance conditions can lead to qualitatively different 
results, similar in stability analyses. For instance, the former condition 
leads to the conclusion that electromagnetic waves may be amplified if the 
energetic electrons possess a loss-cone type distribution and the latter cannot 
achieve this result. This point can be easily envisioned from the analysis 
presented by Wu and Lee (1979). 

We shall return to this point later when we discuss the numeri- 
cal results to be presented in Sec, IV. 

B. Method of Calculating the Dielectric Tensor 

since the dispersion equation can be writter in terms of the dielectric 
tensor € (k, U ) , i. e. 

•I 
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Co obcein che expreesloas for Che censorial coatponeacs of appears 


Co be a crucial Cask. Hereafter we Incroduce Che ieflnitlon 


^Ij “ ^Ij * \j 


(9) 


where che superscripts b and e denoce respectively che concrlbucions 
from che background (or backscatcered) electrons and che energeclc electrons. 
Because Che background electrons are supposed co have low energies, 
may be written as 

2 
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where « 4ne 


These expressions are well known in Che litera- 
ture and obtainable with che so-called "cold-plasma" approxlmaclon. 


On the ocher hand, for the energetic electrons, can be ex- 


pressed as 


OJ f 

Q*j - 2^ *4 I 
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e z z 


T\"^ 's:;' * “u -507) ^ 

2tt ^ I dull jdu^[.i ^ J (u, ^ - u, F,(u) 


2 +« 


0 
.(n) 


X I ^ 

^ yw - tu2 - k„u,i 
a«-» e » ii 


( 11 ) 
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vh«r« n - |•B-./o cl ; w • (4ir« n /m ) 1» ch« plaana frtqutncy of 

th« onorgotic «l«ctrons; F is Chs unpsrCurbod distribution function 

ft 

of thn ensrgotic tlectrons; and the tensor is given by 
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where J ■ • 

In Eq. (11) for convenience we have used the nomentua per unit 

mass u = p/m (m is the rest mass of an electron and p is the 

« C e e -> 

electron momentum) instead of v , the electron velocity. Correspond- 

2 2 k 

ingly, we have Y" (1+u/c) which is the relativistic factor. Since 
we are merely interested in waves with u * , we need only to retain 

n»+l and 0 terms in Eq.(ll). Moreover, because the auroral electrons 
have energies about several keV, it can be shown that the argument of 
Che Bessel function is smaller Chan unity. Thus we may use 

the following approximate expressions! 
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B«cA\ui« tht m««n Altctrotx uiargy la low, cha ralativiatlc tfface la 
txpaccad to ba waak. It la apparant that ooly tha affact appaara In 

tha danomiaator; yw - Q - k„u,.. Thla la important, bacauae the pole 

a " M 


datarmlnas tha raaonanca bacvaan tha wava and alactrona, as dlscuaaad 
aarllar. Physically, this condition controls tha amplification or absor- 
ption of tha radiation which concam i ua. 

To procaad with our dlacxiaalon, wa ramark two points. First, ba- 
cauaa wa ara concamad with tha fast axtraordlnary-moda vavas which 
hava indlcas of refraction smallar chan unity, wa saa chat tha anomalous 
cyclotron raaonanca - k||U^I *0 la Impoaaibla. Thus tha u- 

intagral which appaars in tha exprasaion for , assoclatad with 

Q«>1, 0 , may ba again dlscuaaad In tha "cold-alactron" limit. Sacond, 
Cha tarm aasociatad with n * I can ba studiad by first writing (Tsai 
at al.. 1981) 


yuj 


1 




dc exp (Kyui - - 


k„u„)t] 


Wa find Chat subsaquantly tha u intagratlon can ba parformad, If F 
is a loss-cona distribution as dascrlbad by Eq. (2). Tha details ara 
prasancad in a separata artlcla (Wong at al . . 1981). 
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Sxpreaalon for Q^j ■ D 

After some algebraic manipulations, the aforementioned method 
of calculation leads us to the following expressions for Q*^ (n«l) . 

C'-i’ * ' 7 ^ "p(- r) j 
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(15) 


where 0^ (j • 1,2, .,.4), g^(a -) , 2, ...5) and ?g(3 ■ 1,2, ...6) are 

described in the Appendix). Moreover, we note chat ■ iQ^^ ■ - IQ^^ 

- Q and Q ■ Q ■ IQ ■ - iO 
^xx xz ^zx zy ^yz 

In obtaining these expressions we have introduced the following 
definitions: 
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A. Numerical Method 

Evld«i..ly, ch« disparslon equation is complicated and very difficult 
to solve. In the following we use a numerical method for Investigating 
the solution. Before going further ve point out chat since 


exp 


U(l-IT) 


'^1)^ 1 
p-0 ** r(p+i)(i-iT)’ 


• I C~) 


(17) 


Consequently, the censor Q^^Cn-l) may be expressed in terms of a fimction 


F^(z) which is defined as 


F (z) = - i 


' dT 


(18) 


with q • 0,1,2, As pointed out by Shkarofsky (1966), the function 

F^(z) may be written in terms of the usual plasma dispersion function Z 
Fried and Conte (1961) , according to the relation 


1-1 >1 q- 

-ir^( 2 ) . I (- 1 )“ (-2) ^ X[12‘^(12’5)] for q 2 I 


and 


?^(Z) - ^ Zdz^*) 

Z 


(19) 

( 20 ) 


In deriving the expressions for we have assumed implicitly 

chat Che complex frequency oi has a positive imaginary part. Since 
Che analytic continuation of the Z fxmetion to the lower half of w 
plane is well known, Che validity of the function F^(z) can be extended 
to the entire complex plane. This consideration enables us to solve 
Che dispersion equation by numerical method. (n»-l) may be approxi- 
mated as 
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Furthermore, we need to include Q (n=0) - - u f CJ in the dispersion 

2Z C 

equation, because this term is of the same order of magnitude as the other 
terms. 

B. Results of Computation 

Extensive amount of computation has been carried out. Numerical 
solutions of the dispersion equation are obtained for the four cases men- 
tioned in Sec. II. In the following we present the results ir Figs. 7-12. 

In Fig. 7, the maximum growth rate, u. , is plotted versus 
altitude. It is shown that in case C, the maximum growth rate is highest 
at altitude greater than 5000 km. This is the case in which all low-energy 
electrons are assumed to be removed by the parallel electric field. On the 

other hand, the growth rate for case D is lowest. It is seen that CJ. 

° i.max 

decreases rapidly above 5,000 km. At 8,000 km, it is almost two orders 

of magnitude smaller than that in case C. In Fig. 8, the spatial amplification 

distance, (J. /v (where v is the corresponding group velocity), is 
nmax g g 

plotted versus altitude. The distance is normalized. In case B, the spatial 
e- folding amplification distance is the largest in comparision with that for 
case D, where amplification is least effective. On the other hand, the ampli- 
fication is most effectiv , when case C is considered, as shown in Fig. 9. 
This result is attributed to two basic factors: one is the high temporal 
growth rate and the other is the small group velocity corresponding to the 

^^i, max’ 
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Altitude (km) 


Fig. 7. Normalized Growth Rate Veraua Altitude for Density Models A, 
B, C and D, With and Without Parallel Electric Field. (See 
Tables 1 and 2). 
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Fl^«. F. Normalized Amplification Distance Versus Altitude for Models 
With and Without Parallel P.lectrlc Field 
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Altitude (km) 

Fig. 9. Same as Fig. 8, but Depicting Case C With a Parallel Electric 
Field. 


3i 



Pig. 10. Real and laMginary Frequency Veraua Wave Nuaber for W* 
Nomal Anglea of 6 • 0*, 30*, 45* and 60*, Caae A. 



















Wo now e^c'mlne the details of the solutions for cases A, B, C, and D, 
at 6,000 km where Q is about two hundred kHe. It is found that in 

Q 

case A and case B, the maximum growth rates for various 6 are some- 
what similar, if #£60°. For #>60° the growth rate diminishes. We 
only present the result concerning case A. Both the growth rate and 
the real frequency are plotted in Fig. 10. The upper panel shows 
as a function of k. Four values of S are considered. For each value of 9 , 
the growth rate occurs over a certain range of k. The peak value of 
Uj shifts to higher values of the wave-number as # increses. In the lower 
panel, the real frequency is plotted versus k. 

Case C is presented in Fig. 11. The upper panel displays the growth 

rate and the lower panel shows the real frequency. The abscissa in kc/ Q . 

o 

An Interesting feature is that in this case the cutoff frequency becomes lower 
than the cyclotron frequency. Thus, even if kj| = 0, wave-electron re- 
sonance is still possible. In fact, in this case, the growth rate is higher 
for # = 90°. Except in the region kc<xfl^, the growth rate decreases rapidly 
as # decreases from 90°. 

Lastly, we present case D in Fig. 12. In this case we find that the 
growth rate is significant only for 0 S # < 45°. For S > 45°, the growth 
rate decreases rapidly. The growth rate is computed for three values 
of # , namely ♦ = 0 , # = 30°, and S = 45°. The result is about one order 

of magnitude smaller than that for case C. Furthermore, the bandwidth of 
the unstable mode for each given 0 is fairly narrow. 

C. Discussion 

If we consider a loss -cone as indicated in Fig. 13 in which the 
boundary of the distribution function is shown to intersect with the 
resonan e circle, the necessary condition for this to occur in that 
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Qc = loss -cone angle 


Ftp. 13. Schematic Figure Depleting the Physlca of the Amplification 

Process. The resonance surface Is shown In typical dimension 
with respect to the upcoming atmospheric less cone. Under 
normal circumstances the loss cone region will have negative 
anisotropy and contribute to wave growth. 
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This condicion may ba revrictsn as 




T 2 

cos 6 cos a 


If we fix a and (cc /O - 1) , it is seen Chat as 6 increases 

from 0* , the threshold value of for instability also tends 

to Increase. This explains why the threshold unstable wavenumber tends 

to shift to higher values when 6 changes from 0* to larger values, 

as shown in Fig. 10 and Fig. 12. For case C, the wave frequency may 

be lower than (2 . We see from Eq. (22) that in this case the in- 
e 

stability can occur, even if k • 0 . This observation agrees with 
the numerical result shown in Fig. 11. 

It is true in general, that growth rate diminishes when k becomes 
sufficiently large. This is expected because the larger the k , the 
higher the (o^ - (2^ and, therefore, the fewer the resonant electrons. 
The frequency bands of the unstable modes are fairly narrow for all 
cases studied in this paper, mainly because the characteristic electron 
energies at various altitudes which have been considered in the analysis 
are low. If we have assumed high energies, the growth rates would be 
even higher and the band widths of the unstable modes would be broader 
too. 

The effects of a parallel electric field on the amplification of 
AKR are clearly significant and important. The parallel electric field 
not only modifies the distribution function but also reduces the density 
of the low-energy electrons. Both effects tend to increase the growth 
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rate at a given altitude. However, from cur numerical studies, we find 
that the latter effect (density depletion) prevails over the former. However, 
some time when electron density is fixed, the two distribituon functions 
(1. e. , one with ♦ = 0 and the othev with >1' e 0) can given rise to signi- 
ficantly different growth rate. I’or instance, we have studied density model 
D3 and calculated the growth ratctt for the two cases. The difference be- 
tween the dispersion relations is rather significant. As a result, the growth 
rates have quite different k dcpeovlence. 

The present investigat/ot also finds that the real loss-cone distri- 
bution function is more un'^ ble (i. e, , results in higher growth rates) thf.n 
the iffective loss-sor.e distriuotion function used by Lee et al. , (1980). A 
comparative study 's discussed by Wong et al. . (1981) in which we examine 
the two cases in detail. The effective loss-cone distribution described by 
Eq. (1) peaks in space about Vj^~oi for whatever value of V| , 
whereas the real loss-cone distribution described by Eq. (2) perks in v^ 

space at v., tan a . Thus the populations of resonant electrons for the 
II c 

two cases may be significantly different, if we fix the frequency and the 
wavevector . Consequently, the growth rates obtained are not alike. 

The major conclusions derived from the present study may be sum- 
marised as follows: 

1. ) The real loss-cone distribution function does give growth rates 
higher than that based on the effective loss -cone distribution. 

2. ) Within the context of the physical model and parameter regime 

which have been chosen for the numerical computation, we find that when 
n^ > n^ the growth rates are fairly unlfom for waves propagating within a 
certain cone, say 0 $ 9 ^ ’ where 9 is the angle between the 
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ambient magnetic field, B , and the wavevector, k . For case A or 

case B , 6 is about 60° and for case D (Ei, =0) 0 is about 45°. 

o " o 

On the other hand, if n^ = 0, the growth rate is highest when0»9O°. The 
growth rate diminishes rapidly as 0 decreases. 

3. ) The spatial e- folding amplification distance is found typically 

3 

of several tens of kilometers at 4 x 10 km altitude and increases to a few 

3 

hundred kilometers at 8x10 km altitude for case A and B. For case 
D (no parallel electric field) the amplification distance increases rapidly 
with increasing altitude. 

4. ) If all low-energy electrons are assumed to be removed, the 
amplification can be very effective when a parallel electric field is present. 

3 

We find in case C, the e-folding amplification distance at 7 x 10 Um 
altitude is less than a kilometer. This happens mainly because the local 
group velocity is very small. 

5. ) We have examined case C, but considered the distribution 
function with = 0. In this case we find that the cutoff frequency increases 
to a value very close to the cyclotron frequency and the growth rate de- 
creases significantly from that we have obtained for case C. This implies 
that the effect of the parallel field on the unperturbed distribution is very 
important when the density of the low-energy background electrons is very 
low. 
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